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Abstract 



This paper investigates wave properties of hot plasma in a Veselago 
^ I medium. For the Schwarzschild black hole, the 3 + 1 GRMHD equa- 

tions are re-formulated which are linearly perturbed and then Fourier 
OO ' analyzed for rotating (non-magnetized and magnetized) plasmas. The 

I graphs of wave vector, refractive index and change in refractive are 

used to discuss the wave properties. The results obtained confirm the 
presence of Veselago medium for both rotating (non-magnetized and 
magnetized) plasmas. This work generalized the isothermal plasma 
waves in the Veselago medium to hot plasma case. 

^ I Keywords: Veselago medium; 3 + 1 formalism; GRMHD equations; Isother- 

I mal plasma; Dispersion relations. 
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1 Introduction 

Plasmas are found nearly everywhere in nature. These are electrically con- 
ductive and give a strong respond to electromagnetic fields. Plasma fills the 
interplanetary and interstellar medium and are the components of the stars. 
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To explore the dynamics of magnetized plasma and characteristics of black 
hole gravity when it acts in plasma's magnetic field, the theory of general rel- 
ativistic magnetohydrodynamics (GRMHD) is the most accurate academic 
discipline. Schwarzschild black hole has zero angular momentum so plasma 
present in the magnetosphere moves only in the radial direction. 

The gravity of black hole perturbs the magnetospheric plasma. Rela- 
tivists are always curious to study the effect of these perturbations in the 
black hole regime. Regge and Wheeler concluded that the Schwarzschild 
singularity remains stable when a small non-spherical odd-parity perturba- 
tion is introduced. Zerilli [2] explored the same stability problem by consid- 
ering an even parity perturbation. The behavior of electric field generated 
by a charged particle at rest near the Schwarzschild black hole was discussed 
by Hanni and Ruffini |3]. Sakai and Kawata examined electron-positron 
plasma waves in the frame of two fluid equations for the Schwarzschild black 
hole. Hirotoni and Tomimatsu [5] found that a small perturbation of polo- 
dial magnetic field could highly disturb the plasma accretion. They assumed 
non-stationary and axisymmetric perturbations of MHD accretion onto the 
Schwarzschild magnetosphere. Zenginoglu et al. |S] solved a hyperboloidal 
initial value problem for the Bardeen-Press equation to analyze the effect of 
gravitational perturbation on the Schwarzschild spacetime. 

The 3 + 1 formalism (also called Arnowitt, Deser and Misner (ADM) 
[7]) is much helpful to study gravitational radiations from black hole as well 
as analyzing the gravitational waves. To explore the magnificent aspects of 
general relativity (GR), many authors [8]-[l0] adopted this technique. The 
electromagnetic theory in the black hole regime was developed by Thorne 
and Macdonald [HI [12]. Durrer and Straumann [13] deduced some basic 
results of GR by using this formalism. Holcomb and Tajima [H], Holcomb 
[T^ and Dettmann et al. [TB] studied the wave properties for the Friedmann 
universe. Buzzi et al. [T7] investigated the plasma wave propagation close to 
the Schwarzschild magnetosphere. Zhang [18] composed the laws of perfect 
GRMHD in 3+1 formalism for a general spacetime. The same author [19] also 
described the role of cold plasma perturbation in the vicinity of the Kerr black 
hole. Using this formulation, Sharif and his collaborators [20]- [23] discussed 
plasma (cold, isothermal and hot) wave properties with non-rotating as well 
as rotating backgrounds. 

Metamaterials are artificial materials that have unusual electromagnetic 
properties and Veselago medium or Double negative medium (DNG) is its 
most significant class. This medium has both electric permittivity as well 
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as magnetic permeability less than zero. It is also known as negative refrac- 
tive index medium (NIM) and negative phase velocity medium (NPV). Many 
people |2l]-|28] considered this medium to explore the unusual behavior of 
physical laws and their plausible applications. Ziolkowski |29] studied wave 
propagation in DNG medium both analytically and numerically. Valanju et 
al. [20] found positive and very inhomogeneous wave refraction in this un- 
usual medium. Ramakrishna [21] examined the problem of designing such 
materials which have negative material parameters. He also discussed the 
concept of perfect lens consisting of a slab of negative refractive materi- 
als (NRM). Veselago [32] verified the concepts related to the energy, linear 
momentum and mass transferred by an electromagnetic wave in a negative 
refraction medium. In a recent paper [33], we have discussed the isothermal 
plasma wave properties for the Schwarzschild magnetosphere in a Veselago 
medium. The results verified the presence of this medium for only rotating 
non-magnetized plasma. 

In this paper, we investigate wave properties of hot plasma in the vicinity 
of the Schwarzschild event horizon in a Veselago medium. The format of 
the paper is as follows. In Section 2, the general line element in ADM 
3 + 1 formalism and its modification for the Schwarzschild planar analogue 
is given. Linear perturbation and Fourier analysis of the 3-1-1 GRMHD 
equations for hot plasma is provided in section 3. Sections 4 and 5 give 
the reduced form of the GRMHD equations for rotating (non-magnetized 
and magnetized respectively) plasmas. In the last section, summary of the 
results is given. 

2 3+1 Foliation and Planar Analogue of the 
Schwarzschild Spacetime 

The 3 + 1 split of spacetime is an access to the field equations in which 
four-dimensional spacetime is sliced into three-dimensional spacelike hyper- 
surfaces. In ADM 3 + 1 formalism, the general line element is [IS] 

ds^ = -a^dt^ + -iij{dx' + I3'dt){dx^ + (3^dt), (2.1) 

where the ratio of the fiducial proper time to universal time, i.e., ^ is denoted 
by a (lapse function). When FIDO (fiducial observer) changes his position 
from one hypersurface to another, the shift vector /3* calculates the change 
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of spatial coordinates. The components of three-dimensional hypersurfaces 
are denoted by ■jij = 1,2,3). A natural observer associated with the 
above spacetime is known as FIDO. The mathematical description of the 
Schwarzschild planar analogue is given by 

rfs2 ^ -a^(^z)df + dx^ + dy'^ + dz^ , (2.2) 

where the directions z, x and y are analogous to the Schwarzschild coordi- 
nates r, and 9 respectively. The comparison of Eqs. (12.11) and (12.21) yields 

a = a{z), (3 = 0, 'j^ = 1 {i = j). (2.3) 

3 3+1 GRMHD Equations for the Schwarzschild 
Planar Analogue in a Veselago Medium 

Appendix A contains the 3 + 1 GRMHD equations in a Veselagho medium 
for the plasma existing in the general line element and the Schwarzschild 
planar analogue (Eqs. (l2.ip and (12.21) ). In the vicinity of the Schwarzschild 
magnetosphere, the specific enthalpy for hot plasma is [T^ 

/i = , 3.1 

Po 

where the rest mass density, moving mass density, pressure and specific en- 
thalpy are represented by po, p, p and p respectively. For cold plasma, we 
have p = while for isothermal plasma, p ^ but specific enthalpy is con- 
stant. However, specific enthalpy is variable for the hot plasma. This is 
the most general plasma which reduces to cold and isothermal plasmas with 
some restrictions. This equation shows the exchange of energy between the 
plasma and fluid's magnetic field. It is obvious from the above equation that 
for p to be variable, p must be variable. In this paper, we have used the hot 
plasma along with the Veselago medium in the vicinity of the Schwarzschild 
magnetosphere. The 3 + 1 GRMHD equations (Eqs.(yllO)-(A14)) for hot 
plasma surrounding the Schwarzschild event horizon become 

_ = -V X (aV X B), (3.2) 

V.B = 0, (3.3) 

+ (p + P)f V.^ + (p + p)7V.(aV.V)V 

+ (p + p)V.(aV) = 0, (3.4) 
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3^. + (P + P)7V.V^.-i^^^.}(^^ 



+V.V) + 7V.(V.V)(p + p)- {^6., - ^5.5,) V;,V^ 
= -(p + p)7'a, - p, + -^-(V X B),V.(V X B) - ^{aB)^ 

+ - -^[B X {V X (V X (aV x B))}],, (3.5) 

Ana Ana 

(-1 + ^•^)('^ + P)^' - + 2(P + P)7'(V.a) + ip + p) 
aot a ot 

1 r)R 1 r)R 

^ (V X B).(V X aB) = 0. (3.6) 



Ana 

In rotating background, plasma is assumed to flow in two dimensions, i.e., in 
xz-plane. Thus FIDO's measured magnetic fleld B and velocity V become 

V = Viz)e, + n(z)e„ B = B[X{z)e, + ej, (3.7) 
where B is an arbitrary constant. The quantities A, u and V are related by 

m 

V = + Am, (3.8) 

a 

where is an integration constant. Thus the Lorentz factor 7 = , 
takes the form 

When the plasma flow is perturbed, the flow variables (mass density p, 
pressure p, velocity V and magnetic fleld B) turn out to be 

p = + 6p = + pp, p = p^ + 5p = p^ + pp, 

V = V° + 5V = V° + V, B = B° + 5B = B° + 5b, (3.10) 

where unperturbed quantities are denoted by p°, p, V°, B° while Sp, Sp, 5Y, 
SB represent linearly perturbed quantities. The following dimensionless 
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quantities p, p, v^, fz, K and bz are introduced for the perturbed quan- 
tities 

p = p{t,z), p = p{t,z), V = 5V = v^{t,z)e^ + Vzit,z)ez, 

h = — = b-,{t,z)e^ + bz{t,z)ez. (3.11) 

When we insert these hnear perturbations in the perfect GRMHD equations 
(Eqs.([32D-(l3lD), it follows that 

= -V X (av X B) - V X (aV X 5B), (3.12) 

V.(5B) = 0, (3.13) 

ia(Sp^ + + p)7=V.(i| + V,V)v + (p + p)(V.v) 

a ot aot 

= -2(p + p)72(V.v)(V.V) ln7 - (p + p)j\V .VV) .^r 

+ (p + j9)(v.VlnM), (3.14) 



(p + p)r + ^ J + (P + Ph^v^j - —BrB, 

- B X jV X ^ 

47r^ ^ a dt 



+ -L[B X {V X -^}]. + (p + p)7X,V^- + (p + p) 



x-f%v,,kV^V' + 7Vi(V.V)(5p + M + 7Vi(v.V)(p + p) 

+7\(V.V)(p + p) +7'(2V.v)l^,(V.V)(p + p) - -^{(a55,),, 

Ana 

-ia6B,),,}B' = -{6p)i - 7'{(5p + 6p) + 2{p + p)^\V.Y)}ai 

+ -^{(afi,)j - i(^B,),,}6B^ - {p + p)^\v,V^ + vW^Vk,,V'' 
Ana 

-7'{((5p + Sp)V^ + 2(p + p)^^(y.v)V^ + (p + py}Vij 

-7Vi{(5p + 5p)\/^' + 4(p + p)-f\V.^r)V' + (p + p^jV^^tV'', (3.15) 

f -^^^ + v.V(p + p)y - 1^ + (V.V)(5p + My 
a ot a CM 

+2{p + p)7^(V.V)(V.v) + 2(p + p)7'(v.a) + 4(p + p)7^(V.v) 
(V.a) + 2(5p + 6p)-f\V.Si) + (p + p)7'(V.v) + 2(p + p)-f^ 

(V.v)(V.V) + {6p + SpmV.V) = ^[v.(B.^)V + V.(B.— )v 
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+V.(B.5B)V + V.(5B^)V - v.(B.V)— - V.(B.V)^ 
_V,B.,^|?-V,,B.V,f]-J-,V,B.B,|-V,B.f^,B, 

+— [(v X B + V X 6B).{V X B) + (V X B).(V x 6B)] (3.16) 

4:71 

The component form of these equations, using Eg. (13. Ill) , become 
1 dK 



- uh^,z = {uh^ - Vb^ -v^ + Xvz)V In 



a 



a dt 

-{v,,z - At;,,, - \'v, + V'h, + Vh,^, - u%), (3.17) 

— ^ = 0, 3.18 

a dt ' ^ ^ 

K,z = 0, (3.19) 

155 1 dp . , n^^.ldVx , so 

P—^+P-^ + iP + Pn V{ — ^ + uvx,z) + {p + Pn u 
a ot a at a at 

1 f)ii 

x-^ + iP + P)i^ + l'^>^,^ = -iMp + P) [(1 + 27V^)F' 
a at 

u' 

+27Wu>, + (p + p)[(l -27V)(l + 72u2)- 

XL 

-27^uVy>„ (3.20) 

+ |(p + p)7'(l + 7'^') + ^^uvx^, + {{p + p)^^uV 

^^^.^ ~ "7~(1 + - 1— {"'(1 + + 

47r J 47r 47ra ^ ' 

+-f^u{pp + pp) {(1 + 7V^)V' + j^uVu'} + j^uV{pp + pp 

+p'p + pp) + [{p + p)^\ { (1 + 47V^)W + 4VV'{1 + 7 V^) } 

+ + + 27V2)(p' + p')]v, + [{p + p)^' { (1 + 2f n^) 

(1 + 27V2)\/' - 7W' + 27^(1 + 27^2)^^^'} _ l^(Aa)' 
+72^(1 + 2^\^){p' + p')]v, = 0, (3.21) 



1 dv 
X — - 



7 



X 



■^^o; 2 H — — (1 + '^y^x z + — {SS"-^' ~ + u\{uol 
' 4:71 ' Ana 



An ] a dt \ An 

1 { / , ^ 2/-, . 2 2\ , A^i?^ 

A^n A52,. 52 

+ (1 + 72^2) (p'p + pp') + 72u2(p'p + pP') + [(p + p)74 

X{«2l^'(l + 4721^2) ^ ^y^^^ ^ ^ 272m2))| _ 

47rQ; 

+27Vy(p' + [{p + pW {u'{l + 7V) (1 + 472x^2) 

+2u^\a, + (1 + 272^2) W)} + + 272^(1 

+7V)(p' + p')]^. = 0, (3.22) 

1 o dp 1 o dp 2// /\ 2/~ ~ /- 

~^ ^ a7 + ^ a7 + ^ +p)v^ + u-f {pp^^ + pp^^ + PP 

a ot a at 

+p'p) p— + 2-f^u{pp + pp)a^ + l'^u\pp + pp) + 2(p 

a ot 

+p)7^(m\/' + 2M\/a^ + u'V)v^ + 2(p + p)72(272Mn' + 0^7^ 
+272^2^^)^^ + 2(p + p)7^M^^^a.,2 + (p + p)7'(l + 27'«') 
52 



-[(\/2 + ^2^^xK + (T^' + - Ay (Ay 



+^) ^ _ ^( + X,) ^] - [(y - Xu)v,, + A(kA 



9i ' ' 9i ^ Ana^ 

^2 



+ ^(AA't;, - X'v^ - X'Vb, 
An 

+X'ub^ - Vh,,, + mA6^,^) = 0. (3.23) 

For the purpose of Fourier analysis, the following harmonic spacetime depen- 
dence of perturbation is assumed 

p(t , z) = cie-^("*-'") , p(t, z) = C2e-'("*''^^) , 
v^{t, z) = cge-'^'^*-^'"), ?;^.(t, ^) = c^e-'^-^'-^'\ 
b,{t, z) = C5e-^('^*-'="\ 6^(t, ^) = C6e-'('^*-'="\ (3.24) 

Here k and a; are the ^-component of the wave vector (0, 0, k) and angular 
frequency respectively. Plasma wave properties near the event horizon can 
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be explored by the wave vector which is also used to obtain refractive index. 
We define wave vector and refractive index as follows: 

• Wave Vector: The direction in which a plane wave propagates is 
represented by a wave vector. Its magnitude gives the wave number. 

• Refractive Index: When light travels from one medium to another 
(usually from vacuum) then its ratio between the two mediums is given 
by the refractive index. The change in the refractive index with respect 
to angular frequency decides whether the dispersion will be normal or 
anomalous. 

Using Eq.f l3.24p in Eqs. fl3.17p -f l3.23p . we get their Fourier analyzed form 

Ci{a' + ika) — C3 {(«A)' + ikaX } — c^{aV)' — C6{(«'u)' + lu 
+ikua} = 0, (3.25) 

C5( ) = 0, (3.26) 

a 

C5Lk = 0, (3.27) 

Ci(^p) + C2{^p) + cs{p + p)[^i^u + (1 + -f\^)Lk 
a a a 

-(1 - 2^\'){1 + 7V)- + 27VW] + C4(p + Ph'ii^ 

u a 

+iku)V + u{l + 27V^)\/' + = 0, (3.28) 

ci[p-f'^u{{l + 7^2)^/ ^ ^VW} + 7Vu(p' + Lkp)] 
+C2[rf\{{l + 7V^)\/' + 7VW } + 7Vm(p' + Lkp)] 

+C3[(p + pWiil + 27'm2)(1 + 2-f^V^)V' + (^ + Lku)-f^Vu 

a 

--f^V^V + 27^(1 + 2^\^)uVu'} + 7V(1 + 27 V) (p' + p') 

-|^(A«)' + - ^ku)] + C4[(p + phMil + 47V2) 

xuu' + AVV'il + 7V)} + (p + ph^l + 7V2)(^ + tku) 

a 

Wu{l + 27V2 p' + p') + - ^ku)] 

Ana An a 

B^ 

-cq- — [auu' + a\l + u^) + (1 + u^)Lka] = 0, (3.29) 
Ana 
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ci[pi^{a, + (1 + 7^m2)W + 7Vl/V^'} + 7V(p' + tkp)] 
+C2W{az + (1 + 7V)W + 7='mV\/'} + (1 + ^\^) 

x{p' + ikp)] + C3[(p + p)7^{(l + 7^^i^)(— ^ + ^ku) 

a 

+ 7^u^)(l + 472^2) + 2M72(a^ + (1 + 27^2)^^')} 

+27^^1 + 7V) (p' + p') + ^(Aa)' - - iku)] 

Ana An a 

+C4[(p + p)7^{(^ + iA;M)MV^ + mV'(1 + 47V2) + 2r(a, 
a 

+ (1 + 27V)W)} + 27^mV(p' + p') + ^(— - Lku) 

An a 

J + C6[- { — {Aa) + a A — uA[ua 



47ra Ana 

+ -—{l + u'^)Lk] = Q, (3.30) 
47r 

ci{(— ^7^ + ^^^^7^ + 27x7^0^ + 7\')P + V7^} + C2{( — (1 - 7^) 
+6/cM7^ + + 7^m')p + M7V} + C37^{(p' + p ) + 2 

AB^ 

X (27^mm' + + 272M2a^)(p + p) + {I + 27 V) (p + + 

Ana 

x(A'u - + aA'} + C4[2(p + p)7^{(Mr' + 2uVa^ + mV) + uVik} 

+- — {V - uA)luj - aA'] + cq[- {{V'^ + u^)A + A\/(A\/ + 

Ana Ana 

-aA!u + Lka{V -uA)]=Q. (3.31) 
Dispersion relations will be obtained by using these equations. 



4 Rotating Non- Magnetized Flow with Hot 
Plasma 

In this section, rotating non-magnetized background of plasma flow is as- 
sumed, i.e., B = 0. Thus, the evolution equations (13. 2p and (13. 3p of magnetic 
field are satisfied. Substituting i? = = A and C5 = = Cg in the Fourier 
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analyzed perturbed GRMHD equations ( E qs . f 1 3 . 2 8 1) - f 1 3 . 3 1 1) ) . we have 

ClG^p) + C2(^P) + C3(p + p)[^7'^^ + (1 + l^u'')Lk 

a a a 

-(1 - 27'n2)(l + ^\')- + + c,{p + p)7'[(^ 

M a 

+Lku)V + m(1 + 2-i'^V^)V' + 27WM'] = 0, (4.1) 
Ci[p7^m{(1 + 7V^)V^' + 7Vmm'} + 7Vm(p' + tfcp)] + C2[p7^M 
x{(l + 7V)\/' + 7VW} + 7Vm(p' + + C3[(p + p)7' 

x{(^ + iku)-i^Vu + (1 + 27V)(1 + 27V2)V' - 7W' 

+27^(1 + 2-i^u^)uVu'} + 7V(1 + 2-i\^){p' + p')] + C4[(p + p) 

{7^(1 + 7V^)(^ + Lku) + 7\((1 + 47V^)W + 4\/l^'(l + 
a 

7V2))} + 7'«(1 + 27V^)(p' + p')] = 0, (4.2) 
ci[pi^{a, + (1 + 72^2) W + 72M2y\/'} + -f\^{p' + tkp)] 
+C2[p-f^{a, + (1 + 7^M^)W + 7Wy'} + {p + i/cp) X 

(1 + + C3[(p + phHil + + Lku) + U'{1 + 7V) 

a 

X (1 + 47V) + 2^7^(0^ + (1 + 272^2)^^')} + 27^(1 + 7V) 
X (p' + p')] + C4[(p + p)7'{(^ + + «V'(1 + 47V2) 

+2V{a, + (1 + 272^2)^)} + 27Vr(p' + p')] = 0, (4.3) 
ci{(— + '-^^72 + 2M72a^ + 7V)p + wpV} + C2{( — (1 - 7^) 

+Lku-f^ + 27^0^ + 7^u')p + M7V} + C37^{(p' + p') + 2(27^^' 

+ 2-f^u^a,){p + p) + (1 + 272^2) (p + p)iA;} + Ci{p + p)2-f^ 
{{uV + 2uVa, + u'V) +uVLk} = 0. (4.4) 

4.1 Numerical Solutions 

For the rotating non-magnetized plasma, we use the following assumptions 
to find out the numerical solutions 



C -n .11 /l-(tanh(10z)/10)2 

bpecifac enthalpy: p = y — — - — —. 

Time lapse: a = tanh(10z)/10, 
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• Stationary fluid: 0:7 = 1 with velocity components V = u gives the 
following relation: 07 = 1 ^ 7 = — u"^ — V"^ = 1/a. 

• Velocity components: u = V, x and 2;-components of velocity yield 
u = V = -^. 

• Stiff fluid: p = p = fi/2. 

These assumptions satisfy the GRMHD equations (Eqs. fl3.2p -f l3.6p ) for the re- 
gion 1.4 < z < 10, < a; < 10. A complex dispersion relation [34J is obtained 
by solving the determinant of the coefficients of constants of Eqs.f l4.1() - (l4.4p . 
We obtain a quartic equation in k from the real part of the determinant 

Ai{z)k^ + A2{z, Lo)k^ + As{z, uj)k'^ + A^^z, uj)k + A^{z, w) = (4.5) 

which yields four values of k out of which two are real and two complex 
conjugate. The imaginary part gives a cubic equation in k 

Bi{z)k^ + B2{z, uj)k^ + Bs{z, uj)k + B^{z, w) = (4.6) 

from which one real value of k is obtained and the remaining two are complex 
conjugate of each other. From the real values of k ( (I4.5P and (14. 6p ). wave 
vector, refractive index and its change with respect to angular frequency 
are shown in Figures 1-3. The results deduced from these figures can be 
displayed in the following table. 



Table I. Direction and refractive index of waves 





Direction of Waves 


Refractive Index {n) 


1 


Move towards the event horizon 


n <1 and decreases in the region 
1.4 < ^ < 1.8,0 < w < 10 
with the decrease in z 


2 


Move away from the event horizon 


n <1 and decreases in the region 
1.4 < z < 1.6,0 < < 10 
with the decrease in z 


3 


Move away from the event horizon 


n <1 and decreases in the region 
1.4 < z < 1.7,0 < w < 10 
with the decrease in z 



Figures 1 and 2 indicate normal and anomalous dispersion of waves at 
random points while Figure 3 gives normal dispersion in the whole region 
1.4 < z < 10, < a; < 10. The group and phase velocities are equal in 
magnitudes but opposite in directions in all figures. 
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Surface Plot of of Wave Vector 




Figure 1: Waves are directed towards the event horizon. The dispersion 
found to be normal as well as anomalous randomly. 
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Surface Plot of of Wave Vector 




Surface Plot of Phase Velocity Surface Plot of Group Velocity 




Surface Plot of Refractive Index 



Surface Plot of D[n,fi)] 



-6. 95201 
-6. 95201' 
-6 . 95201' 
-6 . 9520i 




Figure 2: Waves are directed away from the event horizon. Region has normal 
and anomalous dispersion at random points. 
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Surface Plot of of Wave Vector 




Surface Plot of Group Velocity 





Figure 3: Waves move away from the event horizon. The whole region 
normal dispersion of waves. 
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5 Rotating Magnetized Flow with Hot Plasma 

In this general case, plasma is assumed to be rotating and magnetized. 
Fluid's velocity and magnetic field are supposed to lie in x^-plane. The 
respective Fourier analyzed perturbed GRMHD equations, i.e., Eqs.f l3.25l) - 
fl3.3ip are given in Section 3. 

5.1 Numerical Solutions 

We consider the same assumptions for the values of lapse function, velocity, 
pressure, density and specific enthalpy as given in Section 4. The restrictions 
on the magnetic field are as follows 



• FoTU = V and = 0, Eq.([3lD yields A = 1. 

The above restrictions satisfy the perfect GRMHD equations (Eqs. fl3.2p - 
dM])) for the range 1.4 < z < 10, < u < 10. We have Cg = from 
Eqs.(|326])-(|327j). Using above assumptions in Eqs.(|325]) and flfasll -f lOTll . 
we obtain two dispersion relations. The real part gives 

Ai{z)k^ + A2{z, Lo)k^ + As{z, uj)k'^ + A^^z, uj)k + A^{z, w) = (5.1) 

yielding four values of k but all are imaginary. The dispersion relation ob- 
tained from imaginary part is 



which gives two real values of k. Figures 4-8 represent their solutions. 
The following table shows the results obtained from these figures. 




Bi{z)k^ + B2{z, uj)k^ + B^{z, uj)k^ + B^^z, uj)P + B^{z, uj)k 
+Bq{z,uj) = 



(5.2) 
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Table II. Direction and refractive index of waves 





Direction of Waves 


Refractive Index (n) 


4 


Move towards the event horizon 


n < 1 and increases in the region 
1.5 < z < 1.7,0 < a; < 10 
with the decrease in z 


5 


Move away from the event horizon 


n < 1 and decreases in the region 
1A< z < 1.8,0 < cu < 10 
with the decrease in z 


6 


Move towards the event horizon 


n < 1 and increases in the region 
lA<z< 1.7,0 < a; < 10 
with the decrease in z 


7 


Move away from the event horizon 


n < 1 and increases in the region 
1A< z< 1.75,0 < tj < 10 
with the (locToase in 


8 


Move away from the event horizon 


n < 1 and increases in tlic region 
1A< z < 1.8,0 < < 10 
with the decrease in z 



These figures show that group and phase velocities are antiparallel. 



17 



Surface Plot of of Wave Vector 




Figure 4: Waves move towards the event horizon. The dispersion is found to 
be normal. 
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Surface Plot of Wave Vector 






Figure 5: Waves move away from the event horizon. Region has anomalous 
dispersion of waves. 



19 



Surface Plot of wave Vector 




Surface Plot of Refractive Index 




gure 6: Waves move towards the event horizon. The dispersion is normal 
well as anomalous at random points. 
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Surface Plot of wave Vector 




Surface Plot of Refractive Index 




Figure 7: Waves move away from the event horizon. Region has normal and 
anomalous dispersion of waves randomly. 
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Surface Plot of wave Vector 




Surface Plot of Refractive Index 




gure 8: Waves move away from the event horizon. The dispersion is normal 
well as anomalous at random points. 
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Normal and anomalous dispersion of waves can be classified into different 
regions given in the following table. 

Table III. Regions of dispersion 



8 



Normal dispersion 



2 < 2 < 10,2.9 < w < 3 
2 < z < 10,3.2 < a; < 3.6 
2 < ;2 < 10,4.3 < cj < 4.4 
2 < z < 10,5.1 < a; < 5.3 
2 < z < 10,6.5 < a; < 6.7 
2 < 2 < 10,7.1 < w < 7.3 
2 < 2 < 10,8.4 < w < 8.6 
2 < ;g < 10,9.4 < < 9.6 



2 < z < 10, 2 < 6^ < 2.2 
2 < z < 10,3.5 < a; < 3.6 
2 < 2 < 10,4.1 < < 4.3 
1.6 < z < 10,5.1 < w < 5.2 
2 < ;2 < 10, 6 < < 6.2 
2 < z < 10,7.7 < a; < 7.8 
2 < ^ < 10,8.3 < a; < 8.6 
2 < ^ < 10,9.2 < g; < 9.4 



1.5 < ;2 < 10,2 < w < 2.2 
1.5 < ;2 < 10,3.5 < < 3.6 

1.6 < ^ < 10,4 < < 4.2 
1.8 < z < 10,5.6 < a; < 5.7 

2 < ^ < 10,6 < w < 6.2 
2 < 2 < 10,7.7 < w < 7.8 
2 < ^ < 10,8.3 < < 8.6 
2 < z < 10,9.2 < a; < 9.4 



Anomalous dispersion 



2 < ^ < 10,2.5 < w < 2.5 
2 < z < 10,3.7 < w < 3.8 
2 < ;2 < 10,4.1 < cu < 4.2 
2 < z < 10,5.4 < a; < 5.6 
2.1 < z < 10,6.2 < a; < 6.4 
2 < 2 < 10,7.4 < w < 7.6 
2 < 2 < 10,8.1 < w < 8.2 
2 < ;g < 10,9 < cu < 9.2 



2 < z < 10,2.6 < a; < 2.7 
2 < z < 10,3.2 < a; < 3.4 

1.6 < ^ < 10,4.4 < w < 4.5 
1.6 < z < 10,5.3 < w < 5.4 
1.6 < ;2 < 10,6.6 < < 6.8 
2 < z < 10,7.1 < a; < 7.2 
2 < z < 10, 8 < a; < 8.2 
2 < ^ < 10,9.5 < g; < 9.9 



1.5 < ;2 < 10,2.6 <uj<2.t 
1.5<z< 10,3.2 < cj < 3.3 

1.6 < ^ < 10,4.4 < cu < 4.5 
1.8 < z < 10,5.2 < a; < 5.4 
2.1 < z < 10,6.3 < w < 6.4 

2 < 2 < 10,7.1 < w < 7.2 
2 < ^ < 10, 8 < cu < 8.2 
2 < z < 10,9.5 < g; < 9.9 



6 Summary 

In this paper, wave properties of hot plasma in the vicinity of the Schwarzschild 
event horizon for a Veselago medium are analyzed. The 3+1 GRMHD equa- 
tions are reformulated for this unusual medium. Linear perturbation in 3 + 1 
perfect GRMHD equations is considered and their component form is de- 
rived. Dispersion relations arc found by using Fourier analysis technique for 
the rotating non-magnetized and rotating magnetized plasmas. 
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In the rotating non-magnetized background, Figure 1 shows that waves 
move towards the event horizon while Figures 2 and 3 indicate that waves 
are directed away from the event horizon. Dispersion is found to be normal 
and anomalous at random points in Figures 1 and 2 while it is normal in 
most of the region in Figure 3. 

For the rotating magnetized plasma, Figures 4 and 6 indicate that waves 
are directed towards the event horizon while waves move away from the event 
horizon in Figures 5, 7 and 8. Dispersion is normal in Figure 4 while it is 
anomalous in Figure 5. Dispersion is normal as well as anomalous randomly 
in Figures 6, 7 and 8. The value of refractive index is less than 1 and also 
phase and group velocities are antiparallel in all the figures which are the 
significant features of this unusual medium. Thus, the presence of Veselago 
medium is confirmed for both rotating (non-magnetized and magnetized) 
plasmas. 

The comparison of the results for isothermal |33] and hot plasma can be 
summarized in the following table. 



Table IV. Comparison of the results 



Results 


Isothermal Plasma 


Hot Plasma 


Existence of waves 


No waves in rotating 
magnetized plasma 


Waves exist in rotating 
magnetized plasma 


Direction of waves 


Some waves move away 
from the event horizon 


Most of the waves move 
away from the event horizon 


Dispersion 


Normal at random points 


Normal in most of the 
region in Figures 3 and 4 



The difference between our work and previous work is that we have used 
the variable specific enthalpy and previous work has been done using constant 
specific enthalpy. In our work, waves exist in rotating magnetized plasma 
while in previous work, there does not exist any wave in rotating magnetized 
plasma. Also, we have found that most of the waves move away from the 
event horizon while for isothermal case, some waves move away from the event 
horizon. This comparison shows that variation in specific enthalpy effects 
the direction of waves. Dispersion is normal at random points for isothermal 
while for hot plasma it is normal in most of the region in Figures 3 and 4. It 
is interesting to mention here that the properties of a Veselago medium turn 
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out for both rotating (non-magnetized and magnetized) plasmas confirming 
its validity. 

Appendix A 

This Appendix contains the Maxwell equations, the GRMHD equations for 
the general line element and the Schwarzschild planar analogue in a Veselago 
medium (e < 0, /i < 0). In this medium, the Maxwell equations are 




V.B = 0, 



(Al) 
(A2) 



(A3) 



V X B= - 



dE 



(A4) 
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Also, the GRMHD equations take the form 

^ + -(B.V)/3 + ^B = --Vx(q;VxB), (A5) 

V.B = 0, (A6) 

2 



2 4 11 DV^ 

- {£tlij)V^}-p\i + ^(V X B),V.(V X B) 
^:^(«B)5 + ^{<^B,hB^ - ^(B X {V X [V 



Sttq;^ Ana Ana 

X (aV X B) - (B. V)/3] + (V X B). V/3})„ (A8) 

^(mf ) - + 0(mf -P) + :^(m7V^i-^- 

Ut ar 2a 

+Plro)£m + 2/.po7'(V.a) + /xpo7'(V.V) - 

x{^Jip,^'v,v, +irf,j) - ^{y X B).(V X ^) - ^ 

x(V X B).(B X ^) - -^(V X B).V/3 - -^e(V x B) 
or 47rQ; Air 

.(VxB) + -^(VxB).(Vxq;B) = 0. (A9) 

Ana 

For the Schwarzschild planar analogue, /3, 9 and i?t7ij vanish, the perfect 
GRMHD equations become 

— = -Vx(q;VxB), (AIO) 

V.B = 0, (All) 

do (9V 

+ (aV.V)po + Po7'V.— + po7'V.(aV.V)V 

+PoV.(a;V) = 0, (A12) 
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-(|^% - iB,B,)V^„ V' + Po7V.{i| + (V. V)/.} 

= -Po/^tV -Pn+^{y X B)^V.(V X B) - -1^(Q;B)^, 

+ 4^(«^^)'^-^'' - 4^[B X {V X (V X (aV x B))}]„ 

+ V.V) W) - - 1 + 2m7'(V.a) + m7'(V.V) 



(A13) 



(A14) 
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